Beyond mean-field study of excited states: Analysis within the Lipkin model 
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Beyond mean-field methods based on restoration of symmetries and configuration mixing by 
the generator coordinate method (GCM) enable to calculate on the same footing correlations in the 
ground state and the properties of excited states. Excitation energies are often largely overestimated, 
especially in nuclei close to magicity, even when transition probabilities are well-described. We 
analyse here the origin of this failure. The first part of the paper compares realistic projected 
GCM and QRPA calculations for selected Sn isotopes performed with the same effective Skyrme 
interaction. Although it is difhcult to perform RPA and GCM calculations under exactly the same 
conditions, this comparison shows that the projected GCM overestimates the RPA results. In the 
second part of this paper, we compare GCM and RPA in the framework of the exactly solvable 
Lipkin-Meshkov-Glick model. We show that the discretized GCM works quite well and permits to 
obtain nearly exact results with a small number of discretization points. This analysis indicates also 
that to break more symmetries of the nuclear Hamiltonian in the construction of the GCM basis is 
probably the best way to improve the description of excited states by the GCM. 

PACS numbers: 21.10-k, 21.10.Re, 21.60.Jz 
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I. INTRODUCTION 

Self-consistent mean-field methods are one of the stan- 
dard microscopic approaches in nuclear structure theory 
0. At present, they are the only available microscopic 
method that can be systematically applied on a large 
scale for medium and heavy nuclei. The full model space 
of occupied states can be used, which removes any dis- 
tinction between core and valence particles and the need 
for effective charges. This allows the use of a universal 
effective interaction, universal in the sense that it can be 
applied for all nuclei throughout the periodic chart. 

Despite its successes, the self-consistent mean-field 
method has a number of well-known limitations. From 
a conceptual point of view, the mean-field approach is 
designed to describe mainly ground-state properties, and 
gives very limited access to excited states. This is in 
contrast to the microscopic methods that are available 
for light nuclei, like the no-core shell model and the 
shell-model Monte-Carlo, which also describe excitation 
spectra. A systematic way to resolve these problems is 
offered by symmetry restoration and configuration mix- 
ing. Several groups now develop methods going beyond 
a mean-field approach based on the generator coordinate 
method (GCM), either with non-relativistic Skyrme 2, 3] 
or Gogny interactions, or with relativistic Lagrangians 
0, see also and references given therein. The aim 

is to obtain a global description of ground states, includ- 
ing collective correlations which cannot be included in 
a mean-field approach, even at the effective interaction 
level, and of excited states of all nuclei in a single, unified 
method. 

An introduction to the method that we develop along 



these lines can be found in Ref. 0. First applica- 
tions have demonstrated that such a method permits 
to describe the energies of low-lying collective excita- 
tions and electric transition probabilities, in-band and 
out-of-bands. However, in many cases it has also been 
found that whenever states can be grouped into rota- 
tional bands, the spectra obtained with the GCM are 
too spread and that the excitation energies are too high. 
For spherical nuclei, in particular those close to doubly- 
magic ones, the low-energy collective spectra are only 
qualitatively in agreement with the data. This feature 
is well illustrated by a detailed study of the quadrupole 
and octupole modes of nuclei around ^°*Pb in Ref. |^: 
the excitation energies of the first 3~ and 2+ states are 
overestimated by more than 1 MeV. On the other hand, 
the energies of giant resonances were found to be more 
realistic, as was also shown in Ref. Q. 

The origin of this problem is not obvious and it can- 
not be expected to be unique. A possible source of error 
could be the inadequacies of effective interactions to de- 
scribe spectra. It is certainly an appealing feature of 
the projected GCM that the same effective interaction 
can be used to generate the mean-field states and to per- 
form their mixing. However, existing interactions are still 
far from describing all nuclei with a similar high qual- 
ity I^OiQ' Furthermore, they have been adjusted ex- 
clusively at the level of the mean-field approximation and 
on ground state properties, with at most constraints on 
the values of global parameters, like the effective masses 
or the compressibility, deduced from the systematics of 
excited states. There is no guarantee that such effective 
interactions will correctly predict spectra, although, up 
to now, GCM results are encouraging, and, in most cases. 
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in good qualitative agreement with available data. An- 
other source of uncertainty of the GCM conies from the 
choice of the variational space in which the configuration 
mixing is performed. It is usually constructed by intro- 
ducing constraints on one or a few collective variables 
related to the shape of the nuclear density in the mean- 
field equations. Such a choice may be more appropriate 
to describe the properties of ground states than those of 
excited states, for which additional degrees of freedom 
might have to be included. 

Up to now, these problems have not been addressed in 
a systematic way, with the exception of the detailed study 
of ^°^Pb and its isotopes mentioned above but which 
did not include the restoration of rotational symmetry. 
Rather than to repeat the same kind of analysis as in 
Ref. (Si] , we will adopt here another strategy. We will first 
compare results obtained with the same effective interac- 
tion using either the GCM or the random phase approx- 
imation (RPA). The RPA is an alternative microscopic 
method to calculate collective excitations, relying also on 
the existence of a nuclear mean field. The RPA is very 
well adapted to the description of collective excitations 
in spherical nuclei. This will enable us to separate the 
problems related to the effective interactions from those 
related to the method. We will make this comparison for 
two Sn isotopes, the doubly magic ^^^Sn and "'^^'^Sn, per- 
forming the GCM and RPA calculations under conditions 
as close as possible. We will see that the results differ in 
a manner that raises questions about the degrees of free- 
dom that are explicitly included in the GCM. However, 
significant differences between both methods cannot be 
easily eliminated and the comparison between both mod- 
els cannot be fully conclusive. 

One therefore needs a view on the problem from an- 
other perspective. Exactly solvable models constitute a 
very fruitful ground for the test of and comparison be- 
tween many-body methods. They also permit to explore 
new developments at a very limited cost. For this pur- 
pose, we need a model where collective variables simi- 
lar to deformations can be introduced and where a dis- 
cretized version of the GCM can be defined in a way 
similar to that of the realistic applications. The Lipkin- 
Meshkov-Glick (LMG) model, introduced in Ref. Tij| , has 
the required properties. Depending upon the strength of 
the interaction, two different kinds of solutions are ob- 
tained at the Hartree-Fock (HF) level Q: "spherical" 
ones at low values of the strength, and "deformed" ones 
beyond a critical strength. The second part of the paper 
is devoted to a detailed discussion of the LMG model. 



II. COMPARISON BETWEEN THE QRPA AND 
PROJECTED GCM ON SN ISOTOPES 

A. Technical aspects 

Our beyond mean-field method has already been pre- 
sented in details and applied to a large number of nu- 



clei 0,11111111 mill Ell. Let us here summarize those 
of its features which are essential for a critical comparison 
with the RPA. 

The starting point of the method is a set of constrained 
mean-field calculations, with the axial quadrupole mo- 
ment as a constraining operator. The wave functions 
obtained for a discrete set of quadrupole moments are 
projected on good particle numbers and on angular mo- 
mentum. For each value of the angular momentum, the 
projected wave functions are mixed with respect to the 
quadrupole moment by the generator coordinate method 
(GCM), leading for each angular momentum to a collec- 
tive wave function spread over a range of deformations. 

The same effective interactions, the Skyrme SLy4 |2l| 
in the mean-field (particle-hole) channel and a density- 
dependent zero-range force 21] in the pairing (particle- 
particle) channel with a strength of —1250 MeV fm'^ are 
used for the construction of the mean-field wave functions 
and for the calculation of the GCM matrix elements. The 
BCS subspace is limited to an energy range of 5 MeV 
above and below the Fermi level. 

Pairing correlations are a necessary ingredient of a 
GCM calculation: the Hartree-Fock Slater determinants 
corresponding to two deformations, for which the num- 
ber of occupied single-particle states of a given symmetry 
is different, are orthogonal. This feature makes a GCM 
calculation numerically unstable. The problem is cured 
by the partial occupation of single-particle levels due to 
pairing. As BCS pairing correlations collapse whenever 
the density of single-particle levels around the Fermi sur- 
face is low, we use the Lipkin-Nogami (LN) method to 
ensure that pairing correlations are present for all values 
of the quadrupole moment. 



B. Results for Sn isotopes 

Let us first compare GCM and RPA predictions for 
the first excited 2+ state in a doubly magic nucleus, 
^^^Sn, which has already been extensively studied by the 

RPA miiiiii. 

In Fig. ^ the mean-field energy curve is shown to- 
gether with the J = to 4 projected energy curves as 
a function of the quadrupole moment. As usual, the 
projected states are labelled by the deformation of the 
mean-field state which is projected. The topography of 
the projected curves is typical for nuclei which have a 
well defined spherical mean-field ground state . After 
projection, there are two nearly degenerate J = min- 
ima corresponding to two slightly deformed mean-field 
configurations, oblate and prolate. The energy gain is 
moderate, of the order of 2 MeV, with a marginal extra 
gain of 210 keV due to configuration mixing. The collec- 
tive GCM wave functions are plotted in Fig.|21for J = 
and 2. 

In Table m we compare the energies and B{E2) values 
obtained in projected GCM and RPA calculations based 
on the same effective interaction j2^ |25| with the exper- 
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FIG. 1: Projected energy curves as a function of the 
quadrupole moment for ^^'^Sn. The particle- number projected 
mean- field deformation energy is represented by a solid line, 
while dotted, dashed and dash-dotted lines correspond to pro- 
jected mean-field energies for J = 0, 2, and 4, respectively. 



imental data |2a. l27l | . The excitation energy of the 2+ 
state is overestimated by both methods. Within the RPA 
framework, one can show that, for Sn isotopes, the inclu- 
sion of phonon couplings decreases the excitation energy 
of the 2+ state and reduces the B{E2) transition proba- 
bility, although this decrease is less pronounced in ^'^^Sn 
than in lighter isotopes l^^. Note also that the energy 
centroid of the isoscalar giant quadrupole resonance is 
predicted correctly by the GCM at 12.45 MeV, close to 
the value given by the empirical systematics 63 A"^/"^ 
MeV. 

We have performed a similar calculation for ^^°Sn for 
which neutron pairing correlations are present in the 
ground state and for all deformations. Results are com- 
pared in Ta ble llTl with a RPA calculation and with 
experiment [2y|. Again, we find a discrepancy between 
the RPA and GCM predictions for the excitation energy 
of the low-lying 2"*" state, whose energy is significantly 
larger in the GCM calculation than in RPA. The B{E2) 
value obtained within GCM is closer to the experimental 
value than the RPA value. This result is consistent with 
our findings for other systems that transition moments 
are in most cases much better described than excitation 
energies, and suggests that also for nuclei that are spher- 
ical at the mean-field level the geometrical properties of 
excited states are better described by the GCM than their 
excitation energies. This may be related to the fact that 



TABLE L Excitaion energy (in MeV) and B{E2) value (in 
e^ fm*) for up transition to the first 2"*" state in ^^^Sn. 



Method 



Energy (MeV) B(E2T) 



GCM 5.69 630 

Particle-hole RPA 5.13 1370 
Experiment 4.04 1100±300 



TABLE IL Excitation energy (in MeV) and B{E2) value (in 
fm*) for up transition to the first 2"*" state in ^^"Sn 



Method 



Energy (MeV) B(E2t) 



GCM 2.40 1350 

Quasiparticle RPA 1.44 440 

Experiment 1.17 2020±40 



the components of the 2+ states corresponding to two- 
quasi-particle excitations breaking time reversal invari- 
ance are completely missing from the GCM model space 
and that these components have a larger contribution to 
energies than to transition probabilities. 

Since neither RPA nor GCM do satisfactorily repro- 
duce the experimental energies and transition probabili- 
ties for both Sn isotopes, it is tempting to conclude that 
the Skyrme effective interaction is not fully adequate to 
describe low energy excitations in this mass region. 

Besides this problem of interaction, a direct compari- 
son between RPA and GCM results shows a discrepancy 
and a suspicion that the GCM variational subspace as de- 
fined in actual applications is not as rich as the RPA one, 
in particular, at small deformations, where two-quasi- 
particle excitations breaking time-reversal invariance are 
not included in our GCM. Unfortunately, there are sev- 
eral differences between both calculations, which might 
affect their comparison, some of them being hard to elim- 
inate. 

In both cases, the same effective interaction is used in 
the mean-field channel. However, since ^^^Sn is a doubly- 
magic nucleus, there are no pairing correlations in the 
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FIG. 2: Collective GCM wave functions for the ground state 
and the first 2+ state of ^^^Sn. 
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spherical ground state with either a BCS or a Bogoliubov 
treatment of pairing. Therefore, there is no pairing at all 
in an RPA approach, in contrast to the GCM. As soon 
as the quadrupole deformation of the constrained mean- 
field state is sufficiently large (in this case, 500 fm^), 
pairing correlations are present also in the BCS method. 
As already mentioned above, the GCM requires generat- 
ing wave functions which vary smoothly along the collec- 
tive path, which is enforced using of the LN prescription; 
hence, in our GCM pairing correlations are present in all 
states, even in the spherical mean-field configuration of 
i32Sn. 

Another difference between the RPA calculation of 
and our GCM calculation is that the Coulomb and the 
spin-orbit residual interactions are not included in the 
RPA, but are always present in the GCM. The effect of 
these terms has been v ery recently studied by Terasaki et 
al. [2^ and Peru et al. • The inclusion of the Coulomb 
residual interaction raises the 2+ excitation energy by 200 
to 300 keV. 

Some recent studies have also pointed out that the 
RPA has shortcomings which are usually forgotten. As 
shown by Johnson and Stetcu, the QRPA does not re- 
store symmetries exactly |3l). Moreover, in some cases, 
the R PA p redicts poorly the ground state correlation en- 
ergies |32l |. the excitation energies [s^ and B{E2) values 

m 

Still, it seems clear that the GCM overestimates exci- 
tation energies in spherical nuclei close to magicity. On 
the other hand, since it looks difficult to remove all the 
differences between GCM and RPA calculations, the use 
of a model that can be exactly solved seems to be the 
most appropriate way to deepen the present analysis. 

III. EXCITED STATES IN THE 
LIPKIN-MESHKOV-GLICK MODEL 

A. The model 

Lipkin, Meshkov, and Glick introduced an exactly solv- 
able model ^3 1 usually called "Lipkin model" or "LMG 
model" in the literature, that has been widely used to 
test methods of approximation for the nuclear many- 
body problem. 

The model consists of N fermions distributed in two A^- 
fold degenerate shells separated by an energy e. In their 
original paper, two different Hamiltonians were proposed. 
The one which is the most usually studied contains a 
monopole-monopole interaction and is given by: 

H ^eJo~^V{J+J+ + J,J.), (1) 

where V is the interaction strength and Jq, J± are quasi- 
spin operators 0, Q, 01 

1 ^ 



N 

J+ = ^ ] C^pC-p, 

J- = jI (2) 

with the algebra 

[J+,J_]=2Jo, [JoJ±]=±J±. (3) 

The operators c^p and cLj, create a particle in the upper 
or lower shells, respectively, where p labels the N degen- 
erate levels within the shells. The operator Jq measures 
half of the difference between the number of particles in 
the upper and the lower levels. 

The exact wave functions are eigenstates of two op- 
erators, the total quasispin operator — ^{J^J^ + 

J-J+) + Jq with eigenvalue J(J + 1), and a signature 
operator e'^"'", which, for an even number of particles, 
has two eigenvalues equal to ±1. Therefore, as discussed 
in detail in Ref. , the interaction docs not mix states 

which have different eigenvalues of and e*'^'^" and the 
Hamiltonian matrix splits into blocks, which are multi- 
plets in J of order 2J+1. The multiplets separate further 
into blocks of size of J and J + 1 corresponding to the 
two values for the signature. 

To understand the connection between the LMG model 
and realistic nuclear models, it is interesting to identify 
the structure of its eigenstates [s^ . In the limit of van- 
ishing interaction strength = 0, the exact ground state 
corresponds to an independent particle state with all the 
lower single-particle levels occupied, while the exact first 
excited state is given by a Ip-lh excitation on top of the 
ground state, the second excited state by 2p-2h excita- 
tions, etc. 

As the Hamiltonian does not mix states which differ- 
ent J values, the exact wave functions arc linear combi- 
nations of the 2 J -I- 1 eigenfunctions of the operator Jo 
within a multiplet of given J 131 . There is one state of 
each possible np-nh content in each multiplet. The non- 
interacting ground state has Jq = —N/2, and belongs 
to the multiplet with maximum J, i.e. J — N/2. The 
first non-interacting excited state is a Ip-lh state which 
has Jq — —N/2 + 1. Pure Ip-lh states are admixtures 
of states from the J — N/2 and J — N/2 — 1 multiplets, 
while pure 2p-2h states are admixtures of states within 
the J = N/2, J = N/2 - 1 and J = N/2 - 2 multiplets, 
etc for higher np-nh excitations until n = N . 

For small values of the interaction strength V , the 
mixing within the multiplets should be small and low- 
lying levels should have a similar structure as the non- 
interacting ones. By contrast, for large values of V , the 
eigenstates will exhibit a complicate mixing of many p-h 
excitations. Thus, the model exhibits a transition be- 
tween shell-model-like states and collective states. 
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B. Mean-Field Approximation 

In mean-field, or Hartree-Fock (HF), approximation, 
the many-body wave function \a,ip) is given by a Slater 
determinant 



JV 



4pl 



(4) 



characterized by two real degrees of freedom a and (p, 
that will be specified below. The particle- and hole- 
creation operators of the corresponding HF single par- 
ticle basis are given by a unitary transformation among 
the operators corresponding to the non-interacting basis 
EIIM: 




cos(a) — sin(Q;) e '^'^ 
sin(Q;) e^'^ cos(a) 




(5) 



The subscripts and 1 denote hole and particle states, 
respectively. The variables a and ip vary both in the in- 
terval [— 7r/2, 7r/2]. They can be identified as constraints, 
that, due to the low dimensionality of the LMG model, 
map the entire space of mean-field states. New quasi-spin 
operators corresponding to the states with finite a and tp 
can be easily constructed, as described, for example, in 
[3^ . An alternative manner to write the constrained HF 
states will be useful in the context of the GCM. Using 
the Thouless theorem [s^, the normalized constrained 
HF states \a, p) can be obtained from the non- interacting 
ground state, that corresponds to a = 1^ = 0, as 



\a,p) 



(a) exp tan(a) exp(i(/3) 



= 0,</. = 0). 

(6) 

A pointed out by Bhaumik et al. ^\ , the constrained HF 
states of the LMG model can also be formulated in the 
language of coherent states, which allows to make use 
of generating functional techniques to calculate matrix 
elements |3q . 

The constrained HF ground-state energy is a function 
of the variables a and 



where 



eN 



cos(2a) + -xsin^(2Q:) cos(2(p) 



X 



{N -l)V 



(7) 



(8) 



Note that, for a given value of a, the lowest HF state al- 
ways corresponds to ip = Q. The eigenvalues of the single- 
particle Hamiltonian, usually called single-particle ener- 
gies, depend on a only for any mean-field state |a, p). 

One can identify the variable a as a deformation pa- 
rameter. There is a phase transition at x = 1 from a 
spherical [a = <p = 0) to a "deformed" ground state. In 
the latter case, the value of a is obtained by solving the 



equation xcos(2a) = 1. The phase transition and the 
properties of exact and approximated ground states in 
this regime were first discussed by Agassi et al. [s^ . 

While the HF states remain eigenstates of J^, "de- 
formed" HF states break the signature symmetry of the 
exact solutions (which is often called "parity" in the lit- 
erature) for any non-zero value of a. The HF states mix 
the np-nh states with even and odd n within a given J 
multiplet. As a consequence, the constrained HF states 
for non-zero interaction strength contain Op-Oh, Ip-lh, 
2p-2h, 3p-3h etc states. The np-nh components with 
even and odd n can be separated with a projection op- 
erator [3^23l- Due to the simple structure of the LMG 
model with one relevant coordinate only, minimization 
of the energy obtained by projection after variation is 
equivalent to projection before variation |4Clj |. 

The signature symmetry is a discrete symmetry, in con- 
trast to the continuous rotational symmetry broken in 
nuclei with a quadrupole deformation. Li et al. 41] have 
introduced a generalization of the LMG model that can 
been used to test techniques for approximate angular- 
momentum projection, see ^42j and references therein. 
The structure of the original LMG model is closer to 
parity projection in octupole-deformed nuclei j40l |. 

The interpretation of the p> degree of freedom is less in- 
tuitive. It enters the HF states as a phase. It is explored 
by time-dependent HF (TDHF) states, hence a necessary 
ingredient of any dynamical model psL |3| . Using the 
variables a and p, one can form a set of two canonically 
conjugate variables with which, for example, the time- 
dependent HF equations can be transformed to classical 
equations of motion |43 |. 



C. Random phase approximation 

The RPA of the LMG model was formulated for the 
first time in Ref. 13]. The RPA is usually constructed 
on top of the "spherical" HF state ]a = 0, = 0), which 
is the ground state for x < 1. In this regime, the RPA 
phonon creation operator, defined as a superposition of 
all possible Ip-lh excitations, is given by 



1 



(9) 



One assumes that the ground state is the RPA phonon 
vacuum ]0), i.e. Q\0) = 0. The first excited state is given 
by Q^JO) with the normalization condition: 



1. 



(10) 



Profiting from the simplicity of the LMG model, the au- 
thors of 13] have solved Q\0) = exactly. The usual 
way is to solve the RPA equations in the space of Ip-lh 
excitations by linearization. Making use of the equation- 
of-motion approach 45, 46]: 



iE-Eo)mSQ,Q^\0), (11) 
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where E is the absoute energy of the RPA state, and Eq 
the energy of the RPA ground state, one obtains the RPA 
equations: 




{E - Eo) 




(12) 



where A = e and B — ~{N — 1)V = — £X- 

From these equations, the excitation energy of the first 
excited state of the Hamihonian within the RPA is 
found to be 



(13) 



This energy is equal to zero for % = 1, where the system 
undergoes a phase transition, and becomes imaginary for 
even larger values of x- 

The RPA is explicitly constructed as a superposition 
of Ip-lh states; hence, it automatically contains the right 
physics of the lowest excited state in the limit x ^ and 
should be accurate in the limit of small x- 

The RPA correlation energy in the ground state is 
given by: 



eN I, 



(14) 



D. Generator coordinate method 



Continoua GCM 



Most applications of the GCM to the LMG model are 
restricted to a study of the GCM ground state and test 
the correlations in the ground-state In such cases, 
one can take the "deformation" a as a single generator 
coordinate. We are here also interested in the description 
of excited states and will also introduce (/? as a generator 
coordinate. 

One can write the A^-particle GCM wave functions as 
a linear combination of the constrained HF states \oi,ip) 
© with an unknown weight function fk {a, ip) 



Iffc) = / da dififkia,^) |a,^). (15) 



Variation of the energy yields the so-called Hill-Wheeler- 
Grifhn (HWG) equation an integral equation for 



J 



da I dip (^{a',ip'\H\a,Lp) - Ek{a',(p'\a,(p)'^ fk{a,p) 
The two kernels entering Eq. H16(l are the norm kernel 



(^'|a, cos(q!) cos(a') + sin(a) sin(a')e''^'^ '^'^ 



N 



(16) 



(17) 



and the Hamiltonian kernel 

{a',^'\H\a,p>) = 



eN 

— {a\ip'\a,if) 



cos2(a) cos2(a') - s\r?{a) s\r?{a')e'^^^'^-'^'^ + x sin2(a) cos^{a')e'^'v + cos^ia) sm^{a')e-'^^f' 



[cos(a) cos(a') + sin(Q;) sin(Q;') e^f"^""^') 



A set of orthonormal collective wave functions are ob- 
tained by an integral transformation of /j. 

Gk{a',Lp')= I da I dip {a' ,(p'\a,ip)^ fk{a,Lp). 

The HWG integral equation can be solved exactly |l3 |. 
Thanks to the simplicity of the LMG model, the GCM 
with a single generator coordinate a gives already the 
exact solutions of the model. The GCM is build 
on signature-symmetry breaking constrained HF states; 
hence, a priori the GCM wave function mixes Op-Oh, Ip- 



Ih, 2p-2h, 3p-3h etc states. However, one can easily see 
that the signature symmetry is restored by mixing of HF 
wave functions corresponding to ±a with weights equal in 
modulus, which comes out automatically from Eqn. H15|l . 
This situation is similar to realistic applications of the 
GCM when a discrete symmetry, like parity, is broken at 
the mean-field level but not for continuous symmetries 
like rotations. 
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2. Discretized GCM 

In realistic calculations 0, S S 1 the HWG equation 
is solved by discretization of the collective variables 

a ^ 

The discretized GCM equations are obtained by replac- 
ing all integrals in Eqns. H15I18|I by sums over discretiza- 
tion points. The integral equation H16|l becomes a matrix 
equation which can be solved by diagonalization. 

As our aim is to understand why the GCM overesti- 
mates excitation energies in realistic applications, we will 
solve the HWG equation of the LMG model by discretiza- 
tion. Since the continuous GCM permits to find the exact 
eigenstates of the LMG Hamiltonian, the number of dis- 
cretization points must be small enough to avoid a trivial 
reproduction of the exact solution. 

We have solved the LMG model for 30 and 50 particles, 
with a discretization on meshes symmetric around in 
a and and with an odd number of points. The meshes 
in a and Lp have been limited to the region in which the 
collective wave function has a sizable amplitude, typi- 
cally less than half the total range of variation of a. The 
mean-field potentials that are obtained for two values of 
X, above and below the phase transition, are plotted in 
Fig. O The mean-field ground state corresponds to a 
and Lp equal to zero when x is smaller than 1, and to a 
non-zero value of a otherwise. The energy surface is very 
flat for small values of (/?; it is only when a is large that 
the energy increases rapidly with ip. Such a topography 
is representative of the deformation energy curve that is 
obtained for spherical nuclei (x < 1) and for nuclei soft 
as a function of deformation (x > 1). 

In Fig. ^ the exact and GCM energies for the ground 
state and for the first two excited states are plotted as 
a function of the two-body interaction strength x for a 
system of 30 particles. We have used a as the only gen- 
erator coordinate. The GCM equations are solved with 
7 equidistant discretization points chosen in the region 
where the collective wave functions have a sizeable am- 
plitude. In realistic applications of the GCM, the overlap 
kernel is used to define the mesh, the requirement being 
that the overlap between two adjacent points is of the 
order of 0.8. However, in the LMG model, this kernel 
does not depend on the two-body interaction strength. 
Since we want to be as close as possible to realistic ap- 
plications, where the exact solutions are not known, we 
have taken the same mesh for all interaction strengths. 
From the collective wave function obtained for a value 
of X around 0.5, we have chosen points in the interval 
(—0.707, -1-0.707). Since the wave functions of the model 
are either even or odd with respect to a, only four dis- 
cretization points are significant. 

For X < Ij the HF ground-state energy does not de- 
pend on X, as the ground state is always "spherical" . The 
correlations beyond mean-field significantly improve the 
HF result for the ground state and bring it very close to 



the exact value for all values of the interaction strength. 
The situation is less satisfactory for excited states: in 
particular, their excitation energies are far from the exact 
values when the interaction is switched off, while the HF 
method gives the exact energies. The discretized GCM 
becomes more accurate than HF only for x larger than 
0.25 for the second excited state and than 0.45 for the 
first one. 

To understand this surprising result, let us analyse in 
more details how the continuous GCM works in the x = 
limit. Let us first note that, while the exact ground state 
wave function, corresponding to a and ip equal to zero, 
is included in the generating functions \a,(p), the first 
excited states corresponding to pure Ip-lh and 2p-2h ex- 
citations are not. One can easily verify that the weight 
functions that permit to extract the exact eigenstates are 
the Dirac 5{a) distribution for the ground state, its first 
derivative for the first excited state and combination of 
(5(a) and its derivatives for higher excitations. The col- 
lective wave functions, given by equation (|18|l . are regular 
functions but with rather sharp peaks. One of them is 
located at a equal to for even signature states, while 
the odd signature states vanish at 0. Moreover, for very 
small values of x, the collective wave-functions have very 
small amplitudes at the most extreme mesh points, leav- 
ing only a very small number of significant discretization 
points. 

Replacing one of the mesh points by a point close to 
the extrema of the wave functions of the first and second 
excited state for x equal to zero, the discretized GCM 
results become very close to the exact values for the three 
first states and for low interaction strengths. For the 
first excited state, which has a node at the origin, this 
discretization is superior to the one based on the x = 0.5 
ground state wave function up to x = 0.75. 



3. The role of the second generator coordinate 

The influence of a second generator coordinate is 
shown in figure [S] where are plotted the energies of the 
first three states of a system with 50 particles. The cal- 
culations are performed with 7 and 9 points in a, and 1 
or 3 points in p. For the mesh in a we choose equidis- 
tant points in the interval (—0.660, -1-0.660). Since the 
mesh in p> is symmetric and p = Qis always a mesh point, 
there is in practice one active point in ip. All discretiza- 
tions give accurate results for the ground state. With a 
mesh in a only, the energies of the first excited state are 
inaccurate, for all interaction strengths with 7 points, 
and below x equal to 0.4 for a 9 point discretization. 
The second excited state is slightly better described al- 
though the accuracy is still limited for small interaction 
strengths. Adding points in p to the calculation with 
7 points in a corrects the behavior near the origin and 
leads to very accurate results for the ground state and 
the second excited state. The energy of the first excited 
state is also improved although there still remains a small 
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FIG. 3: Mean-field potentials for two values of the interaction strength, x = 0.4 and x = 1-2, for A'^ = 30 particles. 



-^25 



-13.5 



lU 



-15J) 



-16J) 




-2Z0 



-23.5 



-24.0 



-25.5 




FIG. 4: Comparison between the exact and GCM energies 
of the first three states of the LMG model for 30 particles. 
The solid lines are the exact result, while the dashed and 
dotted lines represent the GCM solutions obtained with the 
two different meshes (see text) with seven points in a. The 
HF results are indicated by a dash-dotted line. 



FIG. 5: Comparison between the exact and GCM energies 
of the first three states of the LMG model for a system with 
= 50. Two generator coordinates, a and tp are used in the 
GCM. The exact result are represented by solid lines. The 
dotted and dashed lines lines are GCM results with 7 and 9 
points in q; the dash-dotted lines are results of calculations 
with 7 points in a and 3 points in 



discrepancy for small interaction strengths. The combi- 
nation of nine points in a and three points in leads to 
results indistinguishable from the exact ones. 



4- Comparison between the GCM and the RPA 

The RPA permits to determine excited states but also 
correlation energies in the ground state. The correlations 



that are given by Eqn. 114|l are very accurate and makes 
the RPA ground state energies very close to the exact 
values up to interaction strength equal to 0.9. Beyond 
this value, the GCM with 7 discretization points is more 
accurate than the RPA. 

The GCM and RPA results for the excitation energies 
of the first two excited states are compared to the ex- 
act values in Fig. |H1 It is the accuracy of this energy 
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FIG. 6; Excitation energies of the first two excited states as a 
function of the interaction strength x- Exact results are given 
by solid lines, the RPA results by dotted lines. The dashed 
and dash-dotted lines correspond to GCM calculations with 
the same choices of points as in figure 2] 



difference which is the most interesting in usual appli- 
cations of the GCM. The use of a mesh adjusted for an 
intermediate interaction strength [x = 0-5) and on the 
ground state gives values for the energies of the second 
excited state close to the exact ones when the interaction 
strength is large but it fails for weak interactions. The 
situation is even worse for the first excited state which 
does not have the same symmetry as the ground state. 
In this case, results are close to the exact ones only be- 
yond X equal to 0.75, the error being as large as 60% 
for X close to zero. The results are by far better when 
the mesh is adapted for low-x values, as discussed above. 
Then, even for a small number of mesh points, the exact 
results for both excited states are correctly reproduced 
by a discretized GCM calculation. The RPA results have 
a very different behavior. They reproduce the exact re- 
sults quite well for small interaction strengths. The first 
excited state becomes inaccurate above x equal to 0.7; 
the second one deteriorates more quickly and is worse 
than the bhnd GCM discretization above x around 0.5. 



IV. CONCLUSIONS 



not easy to perform fully equivalent RPA and GCM cal- 
culations, it seems clear that the GCM fails to reproduce 
RPA results for spherical nuclei and that at least part of 
the energy overestimation obtained in previous works is 
due to the way we use the GCM. 

The use of a schematic model does not allow to draw 
final conclusions on the origin of the problems encoun- 
tered in realistic GCM applications, but it can give some 
hints on it. To the best of our knowledge, we have ana- 
lyzed for the first time the application of the discretized 
GCM to the LMG model. To summarize our results, 
one can say first that the discretized version of the GCM 
works remarkably well and permits to reproduce the ex- 
act results with a very limited number of points. Of 
course, the dimension of the LMG model is very limited 
but nearly exact results are obtained for the three first 
states using only approximately 1% of the total num- 
ber of independent vectors of the LMG space. A second 
result is that correlations in the ground state are better 
described than correlations in excited state. We have also 
seen that results are closer for the second excited state 
which has the same symmetry as the ground state than 
for the first one. The fact that an appropriate choice of 
a small number of discretization points permit to obtain 
excellent results, better than the RPA, seems to be an 
artefact of the LMG model, hard to transpose on realis- 
tic cases. On the other hand, it seems encouraging that 
the excited states are well described by the GCM when 
the collectivity due to the generator coordinate is large. 
Finally, the introduction of a second generator coordi- 
nate, conjugate to the first one, also improves the GCM 
results and seem to make them more independent on the 
way the discretization is performed. 

How to transpose the LMG results on realistic cases 
is not trivial. The fact that the LMG mean-field states 
break all symmetries of the interaction certainly rein- 
forces the assumption that to break more symmetries of 
the nuclear Hamiltonian will improve the GCM descrip- 
tion of excited states, in particular by introducing states 
breaking time reversal invariance. This could be done 
in several ways. The most sophisticate one would be 
to introduce collective variables conjugate to the defor- 
mation modes of the nucleus, as in the adiabatic time- 
dependent HF formalism of Villars , Goeke and Rein- 
hard |50| , or of Baranger and Veneroni jsj] . A more eco- 
nomic way to proceed along similar lines is to introduce 
a few states breaking time-reversal invariance which are 
guessed to enlarge strongly the variational space: crank- 
ing constraint or specific two-quasi-particle excitations. 
Work along these lines is in progress. 



In this paper, we have analyzed possible causes of the 
inaccuracy on energies of excited states that are predicted 
in GCM calculations. The study of Sn isotopes has shown 
that part of the overestimation of excitation energies can 
be due to the effective interaction: RPA results obtained 
with the same Skyrme interaction for (^^^Sn) are also in 
disagreement with the experimental data. Although it is 
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